Thickness-dependent spontaneous dewetting morphology of 

^ : ultrathin Ag films 

x> '■ 

1 H. Krishna, 2 R. Sachan, 2 J. Strader, 1 C. Favazza, 3 M. Khenner and 2 ' 4 ' 5 Ramki Kalyanaraman* 



O 

Department of Physics, Washington University in St. Louis, MO 63130 
2 Department of Material Science and Engineering, University of Tennessee, Knoxville, TN 37996 
3 Department of Mathematics, Western Kentucky University, Bowling Green, KY 42101 

v. 

4 Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN 37996 
q ; 3 Sustainable Energy Education and Research Center, University of Tennessee, Knoxville, TN 37996 

> 

m 
m 



(N 
O 

o 



Abstract 



£NJ ■ We show here that the morphological pathway of spontaneous dewetting of ultrathin Ag films on 

o. 

SiC>2 under nanosecond laser melting is found to be film thickness dependent. For films with thickness 
h between 2 < h < 9.5 nm, the morphology during the intermediate stages of dewetting consisted of bi- 

X 

continuous structures. For films 11.5 < h < 20 nm, the intermediate stages consisted of regularly-sized 
holes. Measurement of the characteristic length scales for different stages of dewetting as a function of 
film thickness showed a systematic increase, which is consistent with the spinodal dewetting instability 
over the entire thickness range investigated. This change in morphology with thickness is consistent 
with observations made previously for polymer films [A. Sharma et al, Phys. Rev. Lett, v81, pp3463 
(1998); R. Seemann et al, J. Phys. Cond. Matt., vl3, pp4925, (2001)]. Based on the behavior of free 
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energy curvature that incorporates intermolecular forces, we have estimated the morphological transition 
thickness for the intermolecular forces for Ag on SiO?. The theory predictions agree well with observa- 
tions for Ag. These results show that it is possible to form a variety of complex Ag nanomorphologies 
in a consistent manner, which could be useful in optical applications of Ag surfaces, such as in surface 
enhanced Raman sensing. 

1 Introduction 

Silver (Ag) films and nanostructures have strong plasmonic activity and consequently, are very useful in the 
chemical detection of species via surface enhanced Raman scattering (SERS) ffl|2]. It is known that the 
magnitude of the localized field enhancement leading to increase in Raman scattering is very sensitive to 
the roughness or asymmetry (aspect ratio) of the nanostructures [3]. Therefore, controlling the morphology 
and nanostructure characteristics, of metals like Ag, in a reliable and cost-effective manner is very important 
towards further improving the sensitivity and selectivity of SERS detection. 

One potential approach towards creating complex nanomorphologies in metals is to utilize the sponta- 
neous dewetting of thin films |@]-[9]|. In the classical spinodal dewetting instability, an initially smooth film 
is unstable to height fluctuations because attractive intermolecular forces can exceed the stabilizing effect of 
interfacial tension flUEEE]. ^ s a resu lt, a narrow band of wavelengths can spontaneously grow, eventually 
leading to film rupture and, more importantly, to morphologies with well-defined length scales II 1 211 1 311 . One 
of the important observations in polymer dewetting is the behavior of the dewetting morphology as the film 
progresses from its initially smooth state to a final stable state of particles. It has been observed that below 
a transition thickness hj, the intermediate stage dewetting morphology consists of bicontinuous structures, 
while above, it consists of regularly sized holes 111411151 . This change in morphology has been attributed to 
the form of the intermolecular forces influencing dewetting |[T3l . Specifically, for ultrathin films, with thick- 
ness between 1 and 20 nm, the film thickness-dependent intermolecular forces are made up of a long range 
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attractive component, and a shorter-range repulsive component. The transition thickness can be identified 
from the thermodynamic free energy of the system and, as shown by Sharma and Khanna ifTBl . is located at 
the minimum in the curvature of the free energy. Consequently, the appearance of bicontinuous structures is 
correlated with films whose initial thicknesses lie to the left of the curvature minimum, while the formation 
of holes occur in films with thicknesses to the right of the minimum. In polymers, the the magnitude of hj 
is found to be of the order of a few nm's. More importantly, knowledge of this transition can help guide the 
controlled fabrication of materials with different morphologies and length scales. 

Recently, detailed investigations by various authors have shown that nanosecond laser irradiation of 
ultrathin metal films on non-wetting substrates can initiate a similar dewetting instability, producing ro- 
bust and repeatable patterns with well-defined length scales 1161171 [T614T81 . Previously reported articles also 
showed that metal films can have different morphologies, including bicontinuous structures or holes, and 
in the case of Co on Si02 the transition was found to occur between 3 to 4 nm in film thickness |fT9ll20l . 
In this work, we have explored the dewetting pathway for Ag metal on Si02 substrates in the thickness 
range of 2 < h < 20 nm under pulsed laser melting. Besides the numerous applications associated with 
Ag nanostructures, another important reason to choose Ag for this dewetting morphological study is that 
nanosecond laser heating effects do not introduce any novel dewetting effects, due primarily to the large 
thermal conductivity of Ag |fT8ll . To capture the various morphologies between the flat film and the final 
nanoparticle state, the dewetting was investigated as a function of number of pulses. We observed that Ag 
shows the bicontinuous morphology up to a thickness of ~9.5 nm, while above that it shows holes. Moreover, 
measurements of the characteristic length scale showed the expected h 2 behavior associated with spinodal 
dewetting over the entire thickness range investigated. Importantly, the experimentally observed transition 
thickness hj agrees well with the value predicted by using the sum of possible intermolecular forces for Ag 
on Si02, and is analogous to results for polymer films. The results from dewetting of Ag films show that 
various complex morphologies, potentially useful towards plasmonic and non-linear optical properties, can 
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be robustly fabricated in a repeatable and controllable manner. 



2 Experimental Details 

Ag films with thickness from ~2 to 20 nm were deposited in high vacuum (~ 1 x 10~ 8 Torr) by pulsed laser 
deposition at room temperature onto commercially obtained, optical quality, SiC>2 / Si wafers consisting of 
400 nm thick thermally grown oxide layer on polished 5/(100) wafers ifTTl . The deposition rate was typically 
^ 0.3 nm/min. The energy dispersive X-ray spectrometry (EDS) was used to measure the Ag counts of the 
deposited films in a scanning electron microscope (SEM). The EDS counts were converted into an equivalent 
thickness value by using calibration based on step-height measurements of the film thickness. For every film 
thickness we measured the surface roughness via atomic force microscopy (AFM) and established an upper 
limit of 0.5 ±0.2 nm for the average root mean square (RMS) roughness over the entire thickness range. 
Following the deposition, the films were irradiated in vacuum by a varying number of pulses n from a 266 
nm ultraviolet laser having a pulse length z p of 9 ns. Irradiation was at normal incidence by an unfocused 
laser beam of area lxl mm 2 at a repetition rate of 50 Hz. Under these vacuum deposition and irradiation 
conditions, the film surface was never exposed to air and hence no role of oxygen or an oxide layer was 
expected. The dewetting morphology was investigated as a function of film thickness h, and the number of 
pulses n, which typically ranged between 10 to 10,500 pulses, for irradiation at laser energies at or just above 
the melt threshold of the films. We have shown earlier that the melt threshold energy is a function of the film 
thickness [20]. For each thickness the melt threshold was determined by a visible roughening of the metal 
film surface, as detected under high-resolution SEM within the longest time scale of the experiment (i.e. after 
10, 500 laser pulses) ETTl . The range of E used for the thickness regime investigated here was 60 < E < 120 
mJ/cm 2 . For this irradiation condition, the heating and cooling rate of the metal film was of the order of 10 10 
K/s with a total heating plus cooling time per pulse of ~ 100 ns, which was much smaller than the spacing 
between pulses of 20 ms. Consequently, as we have quantitatively shown earlier [22 ], negligible contribution 
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to the morphology evolution was expected from processes in the solid state. Hence, any morphology changes 
occurred primarily during the liquid phase following each pulse. We also confirmed by EDS measurements 
that the laser irradiation did not result in substantial evaporation of the Ag, even after the longest irradiation 
experiments of 10,500 pulses. 

3 Results and discussion 

Figure [Ja) shows the AFM micrograph of as-deposited 2 nm Ag film on Si02 substrate. The average RMS 
roughness of the film was measured by plotting the profile of the film along a horizontal dashed line (see Fig. 
(Ha)), which is shown in Fig. [2(b), indicating RMS roughness to be ~0.5±0.05 nm. Similar measurements 
were made for all the films investigated. In Fig. [2](a-f), a series of SEM images denoting the early stage 
dewetting morphology, following irradiation by 10 laser pulses, is shown for a different Ag film thickness. 
The important evidence from this series of images is the distinct transition in the general nature of the mor- 
phology between the 9.5 and 11.5 nm films. Up to 9.5 nm [Fig. Oa-d)], the morphology generally consists 
of asymmetric undulations or a bicontinuous type structure. On the other hand, for the 11.5 and 20 nm films 
[Fig. |2£e-f)], regular holes are clearly visible. Quantitative information about the spatial characteristics of 
these morphologies was obtained by evaluating the fast Fourier transform (FFT) of the SEM image contrast, 
which is related to the dewetting film's height variations. The resulting FFT information of the contrast 
correlation is shown in the inset of each figure, and the important information here is the annular form for 
each of the films. This annular FFT is indicative of a narrow band of characteristic length scales for the 
height variations on the surface. This is an important observation given that the dewetting morphology can 
progress via either of the three pathways: homogeneous nucleation, heterogeneous nucleation, or spinodal 
dewetting ||23Tl . In the case of homogeneous nucleation, the features are randomly distributed, both spatially, 
and in time, and no characteristic length scale should appear in this type of dewetting f2M . Heterogeneous 
nucleation can occur due to defects, impurities or other experimentally imposed heterogeneities. In this type 
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of dewetting, a characteristic length scale could appear at the early stages of dewetting only in the pres- 
ence of available ordered nucleation sites. However, we have not observed spatially ordered heterogeneities 
on the substrate surface, as well as on the as-deposited films, prior to irradiation. Therefore, the results 
presented here point strongly to the third option, which is spinodal dewetting ll25l . 

The characteristic length scale associated with spinodal dewetting is established at the very early stages 
of film deformation ffl. The early stage undulations, which occur prior to the appearance of large height 
variations in the film, are extremely difficult to capture experimentally. However, the subsequent morphol- 
ogy, which is a result of ripening of the initial undulations, forms as dewetting progresses and has length 
scales directly related to the initial length scale. Therefore, as shown by many authors, the final nanoparticle 
length scale, can be used as a measure of the thickness-dependent behavior @[TT|]. Here, the length scales 
were measured at different stages of the dewetting process as a function of film thickness. In Fig. [3ja-c), 
the progression of the morphology is shown for the 4.5 nm film, while Fig. [3je-g) shows the same for the 
11.5 nm films as a function of laser pulses between 10 and 10,500 shots. Figures Od) and (h) are the radial 
distribution functions (RDF) for each stage for the 4.5 and 11.5 nm film, respectively. From the position of 
the peaks in such RDF measurements, obtained directly from the FFT's, it was possible to generate the char- 
acteristic length scale present in the pattern at each stage. The result of measuring length scales from these 
progressions is shown in Fig. HJa). The early stage behavior is shown by the closed squares, the intermediate 
stage is shown by open triangles, while the final nanoparticle state is shown by open circles. One important 
observation from this measurement is that no dramatic change in length scale is seen when the morphology 
changes from the bicontinuous to the hole structures, i.e. between 9 to 1 1 nm. This strengthens the argu- 
ment that both morphologies can arise for spinodal dewetting of Ag. Since the characteristic length scale 
for spinodal dewetting is known to vary as X sp °< h 2 , we have also plotted h 2 trend lines for the early stage 
(solid line) and nanoparticle stage (dotted line) data sets with an h 2 trend. The early stage dewetting agrees 
reasonably well with the spinodal trend over the regime investigated, confirming a previously reported result 



6 



for Ag, based on the behavior of the nanoparticle state |[T8l . An important observation can be made on the 
apparent deviation of the intermediate and final stage length scales for the thickest film investigated, i.e. the 
20 nm film. The nanoparticle length appears much smaller than the trend (dotted line), while the interme- 
diate state appears much larger. This can be understood as follows. The drop in the nanoparticle length 
scale is due to the nearest-neighbor interparticle spacing being dominated by a Rayleigh-like break-up of 
the arms of the polygon, as shown in Fig. HJb). In contrast, particles in the 11.5 nm film form at the vertex 
of the polygons formed during the dewetting stage [Fig. [Hg)]. A similar argument, based on a change in 
the feature shape being measured, is likely to explain the intermediate state behavior. As shown in Fig. 
Se-g), the progression of the hole morphology is through merging of the holes into polygons, whose size 
(diameter) will be dependent upon the number of holes it was formed from. It is quite likely that the large 
increase in length scale for the intermediate state of the 20 nm film is because the measured length scale is 
for polygons formed from numerous holes. 

The above results provide a strong case for Ag dewetting via the spinodal process, with additional 
support coming from that fact that similar morphological characteristics are observed when polymer films 
dewet by the spinodal instability. As mentioned in the introduction, polymer films have been observed to 
have a bicontinuous to hole transition at thicknesses hj of a few nm and the theoretical position of this 
thickness has been correlated to the nature of the free energy 11131115] . Specifically, the location of minima in 
the free energy curvature has been correlated to this transition thickness. To determine if a similar behavior 
is seen for the Ag metal, a sum of different types of attractive and repulsive intermolecular interactions have 
been used to estimate the free energy and its curvature as a function of thickness h for Ag on Si02- Since 
the films are very thin, we have neglected the gravitational term in our calculations. 
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3.1 Free energy analysis to determine transition thickness hj 

The intermolecular interaction free energy of a uniform thin film can be realized by first describing the 
disjoining pressure acting on a film. Consider a liquid film of height h on top of a flat substrate. The 
thickness dependent disjoining pressure, 11(A), can be expressed by adding a long-range attractive (Van der 
Waals interaction) and a short range repulsive interaction between the various interfaces formed as a result 
of having the film on a substrate. The repulsive term typically consists of two parts: (i) a Lennard- Jones 
(L-J) type repulsion ll26l which appears due to the electron cloud interaction, and (ii) an electrostatic force, 
which is a result of electric layers forming at the liquid-substrate interface. 

The total disjoining pressure, by considering a long-range attraction, a short range repulsion expressed 
as a Lennard- Jones type form ll26l . and the electrostatic force H3l|26l, is given by: 



n <« = M 



ft,V 1/6, N ' 



s p 

+ — exp{-h/l) (1) 



h J 3\h 

Here, A is the Hamaker coefficient, which has a negative value in units of Joule, and h c is the critical length 
at minimum H(h). S p is the spreading coefficient which is related to the magnitude of electrostatic part 
of the disjoining pressure, and / is a correlation (or Debye) length. Using Eq. [Q the free energy density 
(energy/area) of a uniform film can be written as: 

A Ah^ 

AG= V?-2A^ +SPeXp{ - h/l) (2) 

where, h c is defined as: 

A 3 4 / 3 

-— = -2ysin 2 {e/2) (3) 

and hj (the transition thickness) is calculated from the position of the minimum in the free energy curvature 
as: 



8 



d 2 AG 3A 3Ahf. SP . , n 

^ = F -^ + 7T^(-V0=0 (4) 

The above analysis was performed for Ag on Si02 substrate, and the correlation length / was taken to 
be in the range of 0.2 — 1.0 nm ll27l . The values of the constants used in the analysis for Ag are given in 
Table Q] Table Q] also consists of material parameters for few other metals of interest, where the value of hj 
is estimated by considering that / lies between 0.2 to 1 nm. The typical transition thickness (hj) calculated 
using Eq. @]for Ag appears to be between 0.3 and 11.45 nm. The free energy and curvature plots using 
Eq. |2] and 0] and the correlation length / = 1 nm for Ag are shown in Fig. [5] With this result we can make 
the following observations. In polymer films, the bicontinuous structures are observed for the thinner films, 
while holes are observed in the thicker ones. Similarly, for the presented case of Ag films on Si02 substrate, 
the bicontinuous structures are obtained for films hj < 9.5 nm, while above this thickness the morphology 
evolution starts with the formation of holes. This experimental result matches well with the theoretical 
prediction of 1 1.45 nm, provided a correlation length of / = 1 nm is used. 

4 Conclusion 

In conclusion, we have investigated the morphology evolution in nanometer thick Ag films (2 < h < 20 nm) 
under melting by nanosecond laser pulses. The initial dewetting morphology is found to be a bicontinu- 
ous morphology for films with thickness <9.5 nm, while discrete holes appear for films with thickness > 
11.5 nm. Evaluation of the characteristic dewetting length scales indicated a good agreement with a h 2 trend, 
which was evidence for the spinodal dewetting instability. The observations of the different morphologies 
and a transition thickness is consistent with the behavior observed previously for polymer films. The ex- 
perimentally observed transition value for Ag, of between 9.5 - 11.5 nm, agrees well with predictions from 
intermolecular forces. This work shows that complex and controllable morphologies can be obtained for 
Ag, but more insight is needed to understand the physical origins of dewetting pattern morphologies in thin 



metal films. 
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Table 1 : List of metal parameters required for free energy analysis 



Metal 


YM/Vacc 
(J-m- 2 ) 


e 


SP 
(J-m- 2 ) 


A 
(J) 


h c using Eq. [5] 


hj using Eq. [4] 




m 


(degree) 


na 


El 


(nm) 


(nm) 


Ag 


0.925 


82 


0.79 


-1.33xl0" 19 


0.300 


0.3 - 11.45 


Au 


1.145 


88 


1.11 


-1.76xl0" 19 


0.293 


0.29 - 11.92 


Co 


1.882 


101 


2.26 


-3.1xl0' 19 


0.273 


0.27 - 12.10 


Cu 


1.304 


92 


1.35 


-2.08xl0" 19 


0.288 


0.29 - 11.62 


Fe 


1.870 


101 


2.24 


-2.9xl0' 19 


0.265 


0.26- 11.88 


Mn 


1.152 


88 


1.12 


-8.12xl0" 22 


0.019 


2.64 - 19.23 


Ni 


1.781 


100 


2.10 


-9.9xl0 _2 ° 


0.160 


1.13 - 13.47 


Pt 


1.746 


100 


2.04 


-2.56xl0" 19 


0.259 


0.26-11.94 


Ti 


1.525 


96 


1.69 


-4.32xl0' 19 


0.372 


0.37 - 10.72 


V 


1.770 


100 


2.08 


-1.74X10 -19 


0.212 


0.21 - 12.58 
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Scanning length (urn) 

(a) (b) 

Figure 1: (a) AFM image of an as deposited 2 nm Ag film on SiC>2 (image size is 5x5 fim 2 ). (b) A line 
profile taken along a horizontal dashed line of the AFM image (a) indicating the root mean square roughness 
(RMS) of the Ag film is ~0.5±0.05 nm. 
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(e) (f) 



Figure 2: (a-f) SEM images of the early stage dewetting morphologies following irradiation by 10 laser 
pulses. The film thickness from (a) to (f), corresponds to 2, 4.5, 7.4, 9.5, 11.5, and 20 nm, respectively. 
Also shown in the inset of each figure is the FFT of the contrast correlations of the SEM images. The 
annular spectrum indicates that a well-defined length scale characterizes each pattern. The size of image 
(a) is 0.5x0.5 jKm 2 , images (b-d) are 1.5x1.5 jUm 2 , hid images (e) and (f) are 3x3 jUm 2 and 10x10 jKm 2 , 
respectively. 




Figure 3: (a-c) SEM images (1.5x1.5 fim 2 ) of progression of dewetting in the 4.5 nm thick film with in- 
creasing number of laser pulses (10, 100, 10,500). (d) Plot of the radial distribution function (RDF) for 
each stage of dewetting. The peak position in RDF was used to estimate the characteristic length scales, 
(e-g) Progression (3x3 jl(m 2 ) of dewetting in the 11.5 nm thick film with increasing number of laser pulses 
(10, 100, 10,500). (h) Plot of the radial distribution function for each stage. The letters in plots (d) and (h) 
indicate the RDF's for the corresponding SEM images. 
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(a) (b) 



Figure 4: (a) Plot of the characteristic length scale for various stages of progression as a function of film 
thickness (h). The early stage behavior is shown by the closed squares, the intermediate stage is shown by 
open triangles, and the final nanoparticle state is shown by open circles. Trend lines with h 2 variation for the 
early stage (solid line) and nanoparticle stage (dotted line) data are also shown. For clarity, the error bars 
for the intermediate and final state are not shown, (b) SEM image (15x15 pLm 2 ) of the intermediate stage to 
nanoparticle stage transition for a 20 nm film showing that the break-up of the arms of the polygons is via a 
Rayleigh-type process leading to multiple nanoparticles in each arm. 
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Film thickness (nm) 



(a) 



Figure 5: Plot of the free energy, AG, (solid line) and free energy curvature, ^|r> (dashed line) vs film 
thickness for Ag on Si02 using Eq. [2] and Eq. 0] respectively. The transition thickness hj corresponds to 
the minimum in the curvature and occurs at ~ 1 1.45 nm for a correlation length of 1 nm. 
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